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Abstract. 

Using recent spectroscopic observations, we show that the triple system AO Vel with an eclipsing BpSi primary is 
in fact a remarkable quadruple system formed by two double-lined spectroscopic binaries with components close 
to the ZAMS. All available data have been reanalyzed to derive proper orbital parameters for both binary systems 
and to calculate absolute parameters of the eclipsing system. For the first time, direct determination of the radius 
and the mass have been obtained for a BpSi star. The study of the physical parameters of this unique system is 
especially important since it can be used as a test of evolutionary models for very young stars of intermediate 
mass. 
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fields 


1. Introduction 

Upper main sequence stars with chemical peculiarities 
and large-scale organized magnetic fields constitute about 
10% of A and B-type stars. To understand the physics of 
such stars, it is most important to identify the origin of 
their magnetic field. Two main hypotheses have been pro¬ 
posed: either these stars have retained the primordial field 
present at their formation, the ohmic decay being longer 
than the main sequence stellar lifetime (fossil field theory), 
or the field is generated and maintained by a contempo¬ 
rary dynamo. The first hypothesis implies that strongly 
magnetic stars must be found on the ZAMS (as well as 
over the whole width of the main sequence), while this 
is not necessarily true according to the second hypoth¬ 
esis. Hubrig et al. II2()()()II suggested that magnetic stars 
tend to accumulate in the middle of the main sequence 
strip if their mass is below SMq (which would challenge 
the fossil field theory), while more massive magnetic stars 
seem more evenly distributed over the whole width of the 
main sequence, including near the ZAMS (Hubrig et al. 
Eonu. However, the chemical peculiarities of these stars 

* Based on observations obtained at the European Southern 
Observatory, La Silla, Chile (ESO programmes Nos. 072.D- 
0235, 65.L-0316), and at the Complejo Astronomico El 
Leoncito (CASLEO), Argentina. 


distort their colour indices, so that their effective temper¬ 
atures suffer from an uncertainty which also propagates 
to other stellar parameters like e.g. the surface gravity 
(North I1998II . Therefore, it is most important to deter¬ 
mine the stellar parameters of magnetic stars in the most 
fundamental way, especially as the extent to which their 
evolution within the main sequence is identical to that of 
their non-magnetic siblings is not known. 

The most suitable systems for the accurate assessment 
of fundamental stellar parameters are double-lined eclips¬ 
ing binaries. While a few double-lined spectroscopic binary 
systems (SB2) containing an Ap star of mass below 3 Mq 
are currently known, the rate of close binaries is much 
smaller among magnetic Bp stars (Gerbaldi et al. 119851 
Carrier et al. 120021 North & Debernardi I2004II and only 
two double-lined eclipsing binaries with a Bp component, 
namely AO Vel (= HD 68826) and HD 123335, are known 
to date (e.g., Clausen, Cimenez & van Houten ri995[ from 
here on CCvH; Hensberge et al. I2004II . In this paper, we 
perform a spectroscopic analysis of the eclipsing SB2 sys¬ 
tem AO Vel consisting of two late B-type stars with a BpSi 
primary component (Bidelman & MacConnell ll971?ll . i.e. 
of a type generally displaying rather large spectroscopic 
variations, so that one can infer that there are signifi¬ 
cant surface inhomogeneities. Numerous four-colour pho- 
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tometric observations have been obtained from 1974 to 
1991 with the Danish 50 cm telescope at La Silla (CGvH). 
Times of mid-eclipse for a total of 16 primary and 10 sec¬ 
ondary minima have been determined from all available 
photoelectric observations and photometric elements and 
apsidal motion parameters have been determined from 
complete uvby light curves. This study revealed the pres¬ 
ence of a third component in a wide orbit with an or¬ 
bital period of 25.6 yr. The orbital period of the eclips¬ 
ing system appeared to be 1.58 days and the orbit was 
found only slightly eccentric, with a quite fast apsidal mo¬ 
tion (U=56.8yr). Very recently, Wolf & Zejda (1200511 pub¬ 
lished seven additional times of minima obtained from the 
ASAS-3 database lPoimanski l2002ll and new photoelectric 
UBV observations. As a result, they have improved the 
apsidal motion parameters and the orbit of the third body, 
for which they found a period of 33 yr. 

In spite of the considerable amount of observational 
data collected for decades, absolute dimensions of the com¬ 
ponents remained unknown since no accurate radial veloc¬ 
ity observations had ever been performed for this system. 

2. Observations and spectroscopic analysis 

2.1. Observations 

In March 2004 we obtained five spectra of AO Vel with a 
typical S/N of about 100-150 over three nights with the 
echelle spectrograph FEROS at the 2.2 m telescope at La 
Silla, to obtain radial velocities and, ultimately, the indi¬ 
vidual masses of the components. Their wavelength range 
is 3530—9220 A, and the nominal resolving power is 48 000. 
To our surprise, we found that the spectral line profile in 
the FEROS spectra exhibit a variable and complex struc¬ 
ture indicating four spectral companions at some phases. 
In order to confirm the identification of the four compan¬ 
ions and to derive their respective temporal behavior, we 
obtained six additional spectra in January 2005 with the 
2.1m telescope of the CASLEO and the echelle spectro¬ 
graph REOSC. These spectra cover the range 3600-6000 A 
at a resolution of 13 000 and a typical S/N of about 90. A 
careful inspection of the spectral line morphology, along 
with the photometric ephemeris published by CGvH, al¬ 
lowed us to identify in the spectral line profiles four stellar 
components forming two short-period double-lined spec¬ 
troscopic binaries. In the following, we identify these stars 
as A, B, G and D. The two more massive stars form the 
eclipsing binary (A-l-B) with an orbital period of 1.58 days, 
the primary component being the BpSi star. The other 
spectroscopic binary (G-l-D) has a period of 4.15 days, 
and its components have a similar spectral type B9-A0. 

2.2. Spectra separation and RV measurement 

In order to obtain separate spectra for all four compo¬ 
nents of the system and to measure their radial veloci¬ 
ties (RVs), we applied the iterative method described by 
Gonzalez & Levato cnnai, which was adapted here for 


multiple systems. This algorithm computes the spectra of 
the individual components and the RVs iteratively. In each 
step the computed spectra are used to remove the spectral 
features of all but one component from the observed spec¬ 
tra. The resulting single-lined spectra are used to measure 
the RV of that component and to compute its spectrum 
by combining them appropriately. Starting values for RVs 
were measured in the EEROS spectra by cross-correlation 
against spectra of slowly rotating stars of similar spectral 
type, HD 175640 and HD 196426, which had been observed 
with FEROS in May 2000. The RV of these templates 
were determined by measuring individual spectral lines. 
^From these preliminary RVs we calculated the spectra 
of the four individual components using our high resolu¬ 
tion FEROS spectra. Eor these calculations we considered 
only four of our five EEROS spectra, since in the first 
one, taken at MJD 53070.013, three components appear 
strongly blended, making the determinations of RVs less 
reliable. The GASLEO spectra also were not considered 
because of their lower resolution. During the computation 
of the spectra of the individual components, the RVs in 
those four EEROS spectra were redetermined. 

Once the spectra of the four stellar components were 
obtained, they were used to compute the RVs from all our 
11 observed spectra. Eig.^shows the RVs of the four stars 
as a function of time; the left panel shows the multiple 
morphology of the cross-correlation function. 

Table □ presents the measured RVs for all four stellar 
components. Early observations (before MJD 53373) are 
those of FEROS. The assigned RV errors are formal ones 
calculated within the cross-correlation analysis. The qual¬ 
ity of the RV measurements is also influenced by the uncer¬ 
tainty in the computed spectra, which are used to remove 
the lines of the remaining components from the observed 
spectra. In order to evaluate this contribution to the RVs 
errors, we calculated the RVs using different subsamples of 
observed spectra and in most cases we found that RV dif¬ 
ferences are comparable with their formal errors. However, 
in the spectrum taken on MJD 53070.0187 the RVs of com¬ 
ponents A, B and G are not well defined. The RV compu¬ 
tations converge poorly and depend strongly on the initial 
velocities. In this case we assigned larger errors so that all 
possible solutions are included. 

2.3. Spectral analysis 

To normalize the spectra of individual stellar components 
we scaled the combined continuum by their relative lumi¬ 
nosity. We used the relative fluxes in y and b photometric 
bands provided by GGvH for the star A and the star B, 
and the combined flux of G-l-D. Assuming that the fluxes 
of stars G and D are equal, we obtain relative contribu¬ 
tions to the observed continuum of 0.40, 0.28, 0.16, and 
0.16 for components A, B, G and D, respectively. The un¬ 
certainty of these scaling factors is about 5%, and should 
be considered as a possible systematic error in the inten¬ 
sity of spectral lines. 
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Table 1. Radial velocities for all four companions. 


MJD 

RVa 

kms“^ 

kms“^ 

RVb 

kms“^ 

es 

km 

RVc 

kms“^ 

e-c 

kms“^ 

RVd 

kms“^ 

kms“^ 

53070.0187 

-19 

35 

-7 

20 

-7.3 

15.0 

59.8 

1.0 

53071.0875 

-137 

7 

150 

6 

115.9 

2.3 

-81.9 

1.1 

53071.2551 

-167 

7 

196 

6 

115.3 

2.3 

-83.2 

1.0 

53072.0387 

174 

7 

-164 

6 

54.2 

2.2 

-24.3 

1.1 

53072.2683 

76 

7 

-71 

7 

34.1 

2.1 

9.9 

1.0 

53373.1513 

143 

15 

-157 

6 

15.2 

4.3 

31.2 

2.8 

53373.3296 

108 

8 

-89 

6 

51.2 

3.7 

2.3 

2.9 

53374.3348 

63 

9 

-38 

5 

116.5 

2.9 

-83.5 

3.6 

53374.2101 

-16 

11 

76 

6 

120.0 

3.9 

-80.9 

2.2 

53375.1239 

-42 

15 

47 

6 

49.3 

4.1 

1.8 

3.3 

53375.3616 

-151 

18 

180 

10 

11.5 

2.0 

29.9 

5.9 



Fig. 1. Cross-correlation function (left) and RV curves (right) for FEROS (top) and CASLEO (bottom) spectra. Filled 
triangles, open triangles, filled circles and open circles represent stars A, B, C and D, respectively. 


For each component in the system, we calculated syn¬ 
thetic spectra using the SYNTH code ('Piskunov ll992|l . the 
VALD database (Kupka et. a,1 II 99911 and ATLAS9 model 
atmospheres (Kurucz I1993II . Starting from models with 
different atmosphere parameters, we compared a variety 
of syntheses with observed spectra to determine the best 
fits (TablelSJ. Because of the rather poor quality of the 
spectra, especially of the faint components C and D, we 
mainly took into account the region 4400-4600 A which 


includes the Hei 4471 and Mgii 4481 lines, as well as the 
5000-5100 A and 6300-6400 A regions with some promi¬ 
nent Si- and Fe-lines. The estimated errors in Tes are of 
the order of ± 250 K for component B and ± 500 K for 
components C and D. The best choice of microturbulence 
velocity, was found to be 0 kms“^ for A and B and 
2 kms“^ for C and D. 

Since the deeper eclipse in the primary system, A-l-B, 
corresponds to the occultation of star A, this star is not 
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only the larger and more massive one, but also the hotter 
one. However, He i lines are weaker than in star B, while 
Si II lines are very strong. Since the component A is a 
chemically peculiar star, it is not possible to derive its tem¬ 
perature by simply comparing observations with a variety 
of solar composition models of different T^s and log g val¬ 
ues. A detailed multi-element abundance analysis would 
be needed to determine accurate atmospheric parameters. 
However, such an analysis is not possible yet because of 
the unavailability of high-resolution and high signal-to 
noise spectra. For this reason we adopted for the compo¬ 
nent A the effective temperature predicted by stellar mod¬ 
els for a normal star of the same mass and radius. Using 
the Geneva stellar evolution models (Schaller et al. 1992) 
and the mass and radius calculated from the light and RV 
curves (see next section) we found for the A component 
reff=13700. In Fig.|21we show the observed spectrum (thin 
line) together with the synthetic spectrum computed with 
the atmosphere parameters reff=13700, log 5=4.3 and so¬ 
lar element abundances (solid line). In order to reproduce 
the He i and Mg ii features in the 4470-4480 A region by 
the synthesis using the same atmosphere parameters the 
abundances of He and Mg had to be decreased by 1.45 and 
1.04 dex, respectively (dotted line). From the Sill lines in 
the spectrum of star A we derived a Si-overabundance of 
~ 0.5 dex relative to the solar composition model, which 
is in good agreement with the classification of this star 
as BpSi by Bidelman & MacConnell II1973II . Star B is a 
B8V normal star (Fig.|2). The observed spectrum can be 
very well fitted by a synthetic spectrum assuming solar 
composition, = 12500 ± 250 K and logg = 4.3 ± 0.1. 

The spectra of stars belonging to the spectroscopic bi¬ 
nary C-l-D are similar to each other and correspond ap¬ 
proximately to spectral type B9-A0. The spectral analy¬ 
sis of the system C-l-D using the SYNTH code indicates 
Teff = 10500 ± 500 K for the star C (Fig.EJ, which is the 
more massive companion, and 10000 ± 500 K for the less 
massive companion D (Fig.|2). The Hei line in the spectra 
of C and D cannot be considered as a reliable tempera¬ 
ture indicator since the continuum determination in this 
region is rather uncertain. The uncertainty of our tem¬ 
perature estimates, which are based only on a few well 
resolved, unblended Fe lines, is rather high. The major¬ 
ity of Mg and Fe lines in the spectrum of star C can be 
modelled rather well assuming a purely solar composition, 
reff=10500 K and \ogg = 4.1±0.2. While most metal lines 
in the spectrum of star D can be well reproduced by the 
synthesis assuming solar composition, Teff=10000 K and 
log 5 = 4.2 ± 0.2, the observed Mgii line at 4481 A indi¬ 
cates a lower abundance (w -0.4 dex) than predicted by 
the model. 

For a more careful spectral analysis a larger number 
of observed spectra at different orbital phases is needed, 
not only to improve the signal-to-noise ratio but also to 
assure a better disentangling of spectral components. 


Table 2. Atmospheric parameters derived from compari¬ 
son with synthetic spectra. 


Star 

Tefl 

[K] 

logs 

V sinz 
[kms“'^] 

A 

13700^" 

4.3 

65 

B 

12500 ± 250 K 

4.3 ± 0.1 

73 

C 

10500 ± 500 K 

4.1 ± 0.2 

40 

D 

10000 ± 500 K 

4.2 ± 0.2 

18 


^ Adopted from stellar models according to the observed 
mass and radius. 


3. Orbital analysis 

3.1. Spectroscopic orbits 

Spectroscopic orbits for the two binaries A-l-B and C-l-D 
were calculated using the 11 RV measurements of each 
star. Weights were given to the RV measurements accord¬ 
ing to their uncertainties. 

For the system C-l-D we used the least squares method 
to determine the following parameters: Pen (orbital pe¬ 
riod), Kc, Kd (RV amplitudes), Vycn (center-of-mass 
RV), ecD (eccentricity), ujcd (argument of the periastron), 
and Tcd (time of periastron passage). The RMS of the 
residuals is 4.9 kms“^ for star C, while for star D it is 
1.0 kms“^ for La Silla observations and 4.1 kms“^ for 
CASLEO spectra. These values are consistent with the 
measurement errors. 

In the case of the system A-|-B, we fitted the RV 
curves using the Wilson & Devinney program (Wilson & 
Devinnev ll971[ Wilson lTOOOII in order to take into account 
the proximity effects which are significant in such a short- 
period binary. We lowered the weight of observations taken 
during the eclipses by a factor of 2, since the effect of par¬ 
tial eclipse on the RVs (Rossiter effect, Rossiter 1924) has 
not been specifically modelled during our calculation of 
the component spectra and RV measurements. 

A detailed study of the spectral line profile variability 
during the eclipses would be especially interesting to ob¬ 
tain information about surface inhomogeneities of chemi¬ 
cal elements in the BpSi primary. However, a larger num¬ 
ber of high-resolution spectra are needed for such a study, 
which is thus deferred to a future work. 

Reasonable values for the atmospheric parameters re¬ 
quired by the Wilson & Devinney program were assigned 
according to the spectral types. We adopted the relative 
radii and orbital inclination calculated by CGvH from the 
analysis of the photoelectric uvhy light curves of Grpnbech 
lITWIl . We let the epoch Tab as a free parameter since 
the photometric ephemeris calculated by CGvH were not 
compatible with our RV curves. As described in the next 
section, this and other inconsistencies led us to reanalyze 
the photometric ephemeris and to propose parameters of 
the long-period orbit that are quite different from those 
found by previous works. The orbital parameters resulting 
from the RV curve analysis are listed in TableEl 
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Fig. 2. Observed spectra of components A, B, C, and D (thin lines) compared with synthetic spectra (solid and dotted 
lines). Atmospheric parameters for modeled spectra are as following. A) Solid line: synthetic spectrum with parameters 
from stellar models (Teff=13700 K, logg=4.3) and solar composition, dotted line: same model but with Mg and He 
abundances decreased by 1.45 and 1.04 dex; B) Solid line: synthetic spectrum with parameters from stellar models 
(Teff=13200 K, log(?=4.3), dotted line: best fit for Mg/He ratio and metal lines (Teff=12500 K, log(7=4.3); C) Solid 
line: Teff,min=10000 K, log(7=4.1, dotted line: reff,max=11000 K, logg=4.1; D) Solid line: reff,min=9500 K, logg=4.2, 
dotted line: Teff_„iax=10500 K, \ogg=A.2. 


3.2. The long-period orbit between the two binaries 

To better constrain the fundamental parameters of the 
AO Vel system, we did an in-depth revision of all published 
data and analyzed all available times of minima of the 
eclipsing binary. In the study of the eclipsing system A-l-B, 
CGvH fitted the epochs of 26 eclipses measured through¬ 
out 26 years using a period of 1.584653 ± 0.000005 days 
and considering the combined effect of the apsidal motion 
of the orbit at a rate of 0.0275°/cycle and the light time ef¬ 
fect due to the orbit with a third body (our binary system 
C-l-D). However, the time of minimum predicted for the 
epoch of our observations differs from our spectroscopic 
ephemeris by 0.060 ± 0.026 days. In addition, according 
to the “third body” orbit found by CGvH for the time of 
our FEROS observations the barycentric RV of the binary 
A-l-B should be higher than that of the binary G-l-D for 
La Silla and GASLEO observations, which is absolutely 
incompatible with the observed RVs since, according to 
our measurements, V^ab — Vjcd = —12.3 ± 2.4 kms“^. 
Finally, the third body orbit of GGvH is unable to match 
the 38 old photographic observations taken by Oosterhoff 
& van Houten between the years 1902 and 1935. 

GGvH have admitted this disagreement in their paper and 
speculated that might be a consequence of some kind of 
longer time scale period variation. 


On the other hand, the parameters recently published 
by Wolf & Zejda (2005) for the wide orbit are very different 
from those of GGvH. The ephemeris calculated by Wolf 
& Zejda (2005) agree with our observations within 0.003 
days. They found an eccentric orbit (e = 0.21) with a 
period of 33.3 yr, which might account for the observed 
RV difference between the two binaries. However, it is also 
in disagreement with the old photographic data. 

For these reasons, we decided to recompute the apsi¬ 
dal motion and the third body orbit using all the available 
times of minimum along with our spectroscopic observa¬ 
tions. In view of the very different quality of the various 
photometric sources, we have given weight 1 to the pho¬ 
toelectric minima (26 from GGvH and one from Wolf & 
Zejda 2005), weight 0.1 to the ASAS-3 data, and weight 
0.01 to the much less precise photographic data. We have 
also included the time of minimum predicted by our spec¬ 
troscopic orbit as a single point with unit weight. 

The key difference with respect to the GGvH calcu¬ 
lations is that we do not assume the wide orbit to be 
circular. On the other hand, we have included the photo¬ 
graphic data, which were not considered by Wolf & Zejda 
(2005). Gonsequently, we were able to find a solution that 
is compatible with both our spectroscopic orbit and the 
old photographic observations. Figure 6 shows our fit of 
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Table 3. Spectroscopic orbits of binary systems C+D and 
A+B. 


Parameter 

Units 

Value 


PcD 

days 

4.15008 

± 

0.00016 

Tcd (periastron) 

MJD 

53074.059 

± 

0.060 

V'ycD 

kms“* 

21.2 

± 

0.7 

Kc 

kms“* 

98.3 

± 

2.1 

Kd 

kms“* 

107.4 

± 

1.3 

OJCD 

deg 

60 

± 

6 

ecD 


0.047 

± 

0.015 

acD sin ico 

R© 

16.84 

± 

0.23 

qcD 


0.916 

± 

0.024 

Me sin® icD 

Mq 

1.94 

± 

0.07 

Md sin® icD 

Mq 

1.77 

± 

0.08 

Pab 

days 

1.584584“" 



Tab (periastron) 

MJD 

53069.8821 

± 

0.0055 

Tab (Min I) 

MJD 

53070.8440 

± 

0.0055 

V~fAB 

kms“* 

8.9 

± 

2.3 

Qab 


0.931 

± 

0.033 

UJAB 

deg 

236.5“" 



to 

deg yr“* 

0.0174“" 



eAB 


0.0741“" 



iAB 

deg 

88.5'" 



Ma 

Mo 

3.63 

± 

0.18 

Mb 

Mo 

3.38 

± 

0.18 

&AB 

R© 

10.94 

± 

0.19 

Ra 

R© 

2.34““ 

± 

0.08 

Rs 

R© 

2.1V 

± 

0.08 


Calculated from photometric ephemeris analysis fSec. l8.2(l . 
Pab and uiab are calculated for MJD 53070. 

^ Adopted from CGvH. 

Computed using the photometric relative radii published 
by CGvH. 


the times of minima variation due to the 3rd body (cf. 
with Fig. 4 in CGvH and Fig. 6 in Wolf & Zejda 2005). 
The RMS of the residuals was 4x 10“'* days for the photo¬ 
electric minima, 7 x 10“^ days for the ASAS-3 photometry 
and 2.6x10“^ days for the photographic data. 

In addition to the determination of the parameters of 
the orbit AB-I-CD, the ephemeris fitting provided simulta¬ 
neously the orbital period of the eclipsing binary and its 
apsidal motion rate. Table^l shows the results obtained. 
The parameters we found for the apsidal motion are in¬ 
distinguishable from those found by CGvH and Wolf & 
Zejda (2005). However, the parameters defining the or¬ 
bit AB-I-CD are essentially different from those found by 
CGvH (Po = 25.6 ± 2.5 yr; aopsinf^ = 0.041 ± 0.006 
days) or by Wolf & Zejda (2005) (Pq = 33.3 ± 0.5 yr; 
Oopsinio = 0.052 ± 0.004 days), and also the binary pe¬ 
riod Pab differs significantly. 

With the new parameters for the wide orbit the eclips¬ 
ing binary is at present moving toward the observer, as 
expected according to the observed RVs. The predicted 
RV for the eclipsing binary is —6.0 kms“* with respect 
to the center-of-mass of the quadruple system, which 
is consistent with the observed velocities (see Sec. Oil . 



Fig. 3. Times of minimum shifts caused by the third body 
orbit. Symbols are the differences between the observed 
times of minima and the times of minima calculated using 
mean ephemeris, after correction for apsidal motion (pri¬ 
mary eclipses are marked by squares, secondary eclipses 
are marked by crosses and the filled circle shows Min I 
calculated by fitting our RV curves). 

Table 4. Orbital parameters derived from the analysis of 
times of minima. 


Parameter 

Units 

Value 


Orbit A-l-B 

Pab 

days 

1.5846154 

± 

0.000002 

Tab (Mini) 

MJD 

45043.6608 

± 

0.0004 

P&ps. motion 

yr 

54.72 

± 

0.45 

Tb = 0 

MJD 

39490 

± 

90 

eAB 


0.0741 

± 

0.0004 

Orbit AB -b CD 

Po 

yr 

41.0 

± 

0.2 

To (periastron) 

MJD 

45839 

± 

11 

To (conjunction) 

MJD 

47803 

± 

120 

ao,i sin io 

days 

0.0642 

± 

0.0004 

ao,i sin io 

AU 

11.15 

± 

0.07 

Go 


0.291 

± 

0.005 

iOo 

deg 

11.8 

± 

1.4 

However, more 

accurate 

RVs (errors 

of 

about 0.3- 


kms“* in VyAB — Vjcd) would be required for the mass- 
ratio Mc+D /Ma+b to be calculated with reasonable pre¬ 
cision (5-10 %). 

4. Physical parameters 

4.1. Stellar parameters and evolutionary state 

For the eclipsing pair we have derived absolute masses 
and dimensions using the photometric elements (inclina¬ 
tion and relative radii) published by CGvH along with our 
spectroscopic data. In the mass-radius diagram shown in 
Fig.Elboth A and B are located near the ZAMS. For com¬ 
parison, isochrones of the Geneva stellar models (Schaller 
et af. n 992|l for solar composition are also plotted in Fig.01 
The mass and radius of star A indicates that the age of 
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the system is less than about 5 x 10^ yr, being compatible 
with the ZAMS. This corresponds to less than 25 % of the 
main sequence lifetime for star A. 

Since the orbital inclination for the binary C+D is un¬ 
known, the individual masses of these stars cannot be de¬ 
rived from the RV curves. However, some information can 
be obtained from the outer orbit between the two systems 
and the mass of the binary A-l-B. In the following we use 
subindex “o” for the parameters of the orbit binding both 
binary systems. In particular. Mo is the total mass of the 
quadruple system and qo = yic+D/^ a+b- 

,^From the ephemeris analysis of the eclipsing binary 
we have found the parameter Oopsinio, which is related 
to the semiaxis of the relative orbit (uo) through: 

Oosinzo = Oo,! sinzo • (1 + qo^)- 
Combining this expression with the third Kepler equation 
Mo = Ma+b ■ (1 + Qo) = o^o ■ Pq 

we obtain: 

= (Oo,! sinzo)”^ • Po ■ Ma+b ■ sin^ Zo, 

where (oopsinzo) is given in AU and Po in years. For a 
given value of io this equation provides qo, and hence the 
absolute masses of the system C-l-D: 

Md = qo ■ MA+B/{qcD + 1), Me = qcD ■ Md- iFrom 
the condition io < 90° a lower limit for the masses of 
components C and D can be established. The resulting 
minimum masses are 2.55±O.O9M0 and 2.34±0.10Mq. 

Using these minimum masses for the stars C and D 
and their spectroscopic orbit, we obtain the orbital incli¬ 
nation of the non-eclipsing binary to be icD < 66.0 ± 2.4 
deg. If the mass of C-l-D is close to their lower limit, then 
qo ~ 0.70, and the expected RV difference between the 
two binaries is VpAB — V^cd = —14.5 ± 0.6 kms“^, in 
agreement with the observed value (—12.3 ±2.4 kms“^). 

The parameters of the wide orbit restrict the position 
of stars C and D in the mass-radius diagram to the right 
of the two continuous vertical lines in Fig.^l 

Assuming the four stars to be coeval, stars C and D 
would be very close to the ZAMS, considering the youth 
of the more massive stars. Even the possibility of them be¬ 
ing pre-main sequence stars must be considered. In Fig.^] 
pre-main sequence isochrones for solar composition from 
Bernasconi iHTOll are plotted. The spectroscopic tempera¬ 
tures and the luminosity difference between the two binary 
systems provide useful constraints on the evolutionary 
state of the stars C and D. Three isotherms have been in¬ 
terpolated in the pre-main sequence models of Bernasconi 
Ill99till at temperatures 9500, 10000, 10500, and 11000 K. 
These curves define in Fig.^Jtwo overlapping strips cor¬ 
responding to stellar models with temperatures close to 
those inferred from the spectra of components C and D 
(9500-10500 K and 10000-11000 K, respectively). Within 
these strips, the younger massive models are not compati¬ 
ble with the observations for various reasons. First of all, if 



Fig. 4. Mass-radius diagram. Filled circles with error 
boxes are stars A and B belonging to the eclipsing sys¬ 
tem. Minimum masses for stars C and D are marked with 
vertical bars. Main sequence (Schaller et al. I1992L con¬ 
tinuous lines) and pre-main sequence (Bernasconi 119901 
dotted lines) isochrones are plotted and labeled with 
log(age). Thick lines are isotherms interpolated for pre- 
main-sequence models. 


M(C) 



5.6 6.0 6.4 6.8 

log age 

Fig. 5. Expected temperature difference of stars C and 
D as a function of age, according to pre-main sequence 
models of Bernasconi imm . 

all four stars are assumed to be coeval, pre-main sequence 
ages of the order or below 10®'^ yr should be discarded 
because of the small radius of star B. Second, massive 
models for stars C and D would be too luminous with 
respect to components A and B. More precisely, models 
with log(age)<6.55 are not compatible with the light-ratio 
derived from the light curve analysis by CGvH. Finally, 
in the upper portion of the temperature strips in Fig. ^ 
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two stars with a fixed mass-ratio Mo/Me = 0.916 would 
have a temperature difference Tc — To larger than ex¬ 
pected. This is a consequence of the rapid evolution at 
larger masses and holds also for the difference in luminos¬ 
ity or radius. As illustration, Fig.^shows the temperature 
difference Tq — To as a function of age. Each curve corre¬ 
sponds to the temperature difference of two pre-main se¬ 
quence stellar models with a given mass ratio (as labeled) 
and with an average temperature {Tc + To)/2 = 10250 K. 
Three mass-ratio values have been considered to take into 
account the observational error of qcD- The similarity of 
the spectral types of stars C and D favors again the older 
models. In conclusion, stars C and D are most probably 
on the ZAMS or only slightly younger. The age of the sys¬ 
tem is not younger than about 10® ® yr in the pre-main 
sequence grid of Bernasconi II1996II . 

Additional spectroscopic observations would help to 
determine the age with higher precision by improving the 
calculation of stellar parameters for the eclipsing system 
and the temperature estimates for the less massive binary. 
In particular, if the pre-main sequence nature of stars C 
and D was established, then the post-main sequence evo¬ 
lutionary age of the peculiar star A would be less than 
about 4% of its main-sequence lifetime. 

An additional argument in favor of the youth of the 
AO Vel system arises from the apsidal motion rate. This 
rate can be used in combination with the stellar and or¬ 
bital parameters to derive the average value of the internal 
structure constant k 2 of the components. The influence 
of the system C-l-D on the apsidal motion of the binary 
A-l-B is negligible due to the long period of the outer orbit. 
Our calculations give logfc 2 = —2.25 ±0.05. This is essen¬ 
tially the same value found by CGvH, since this parame¬ 
ter depends strongly on relative radii but only slightly on 
the masses. Theoretical models for the solar composition 
(Claret & Cimmez ll992|l predict logfc 2 = —2.18 ±0.02 for 
ZAMS stars of the same mass. 

4.2. Absolute magnitude and distance 

Even though absolute masses and dimensions are not pre¬ 
cisely known for stars C and D, their contribution to the 
total luminosity of the quadruple system is rather small, so 
that the absolute magnitude of the latter can be estimated 
with reasonable precision. We have derived the absolute 
visual magnitudes of the components by interpolation in 
the stellar model grids according to their masses. If all four 
stars are assumed to be close to the ZAMS the integrated 
absolute magnitude of the quadruple is = —0.82±0.12. 
This magnitude would be somewhat brighter if the sys¬ 
tem was more evolved (M„ = —0.90 for logr = 7.5) or 
if components C and D were assumed to be still PMS 
stars {My = —0.95 for log tpms = 6.6). Using the appar¬ 
ent magnitude V = 9.34, we obtain an apparent distance 
modulus of 10.22 ± 0.16. 

AO Vel is located in the region of the Milky Way popu¬ 
lated by several OB associations and sparse open clusters: 


open cluster Or 173, OB association identified by Brandt 
et al. iiwni and later referred to as Vel OB2, the “Vela 
sheet” (Eggen 119801119831 Loktin et al. 119971 Loktin & 
Beshenov Ennj behind which lie several young associa¬ 
tions embedded in a dust cloud (Kaltcheva & Hilditch 
2000). AO Vel might be a member of the Vel OB4 associa¬ 
tion identified by the latter authors (and probably related 
with the Vel OB2 association of Slawson & Reed 1988), 
the mean distance of which is 0.72 kpc. 

We note that if the distance to AO Vel is of the order 
of 700 pc, this system might be resolved astrometrically 
as a double star. The maximum separation between the 
two spectroscopic binaries would be then about 27 AU, 
corresponding to about 0.04arcseconds. This is too small 
for the Hipparcos mission — which took place at the time 
of conjunction (MJD 47700), not of quadrature — but 
accessible to interferometers like the VLTI. The maximum 
projected separation occured on MJD 52740 (April 2003), 
so that further astrometric observations of this multiple 
are strongly encouraged. 

5. Conclusions 

Using our recent spectroscopic observations we discovered 
that the triple system AO Vel with an eclipsing BpSi pri¬ 
mary is in fact a remarkable quadruple system close to 
the ZAMS. The analysis of our data as well as all avail¬ 
able published data revealed completely different orbital 
parameters for this system compared to previous works. 
We found that both spectroscopic binaries move in a wide 
eccentric orbit (e = 0.29) with a period of 41 yr and this 
orbit is fully consistent with the observed difference be¬ 
tween the center-of-mass radial velocities of the two spec¬ 
troscopic binaries. The calculated orbit imposes a lower 
limit to the total mass of the binary C±D: McTM^i > 
4.9 Mq. The analysis of the stellar parameters and of the 
evolutionary state of the components reveals their extreme 
youth. Since no magnetic held measurements have been 
performed for this system, we have applied for VLT time 
to obtain spectra in circular polarized light with FORS 1, 
which is equipped with polarization analyzing optics. 

We also plan to obtain additional high resolution, high 
signal-to-noise spectra of this system at different orbital 
phases, in order to improve the determination of the stellar 
parameters and atmospheric elemental abundances of all 
companions. This unique system is especially important to 
test evolutionary models of very young stars of interme¬ 
diate mass. Detailed observations during primary eclipses 
will provide strong constraints on the distribution of the 
abundance anomalies on the surface of the A component, 
complementing the standard Doppler imaging technique 
based on rotation. 
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